Magnetoplasmons and SU(4) symmetry in graphene 
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Abstract. We study magnetoplasmons or neutral collective excitations of graphene in a strong 
perpendicular magnetic field, which can be modelled as bound electron-hole pairs. The SU(4) 
symmetry of graphene arising from spin and valley pseudospin degrees of freedom is explored 
r/2 using Young diagrams to correctly predict the degeneracies of these excitations. The multiplet 

D structure of the states is identical to that of mesons composed of first and second generation 

£h quarks. 

e 

1. Introduction 

Graphene is a single layer of graphite or set of carbon atoms sp 2 bonded together to form a 
honeycomb lattice. The electron in the remaining p z orbital is delocalised and responsible for 
graphene being an excellent electrical conductor [I]. Graphene is associated with numerous 
other superlatives such as strongest material [2] and of course thinnest material; it is the first 
example of an atomically thick membrane. Such properties make it an excellent candidate for 
industrial applications, although in many cases there are obstacles to be overcome before this can 
be put into practice (3]. Aside from its potential usefulness, graphene is scientifically interesting 

{Nl in its own right. Unlike other materials, including the 2D electron gas (2DEG), electrons in 

graphene have a linear dispersion relation close to the Fermi energy and in this sense behave 
quasirelativistically [4j . Effects such as Klein tunnelling, which were previously the dominion of 
high-energy particle physics, have now been observed in graphene [5]. For an undoped system, 
the Fermi energy is located at zero energy. The Fermi surface consists of six points which 
coincide with the six corners of the hexagonal Brillouin zone. Only two of these points are 

^ inequivalent and are commonly termed the K and K' valleys. They may also be referred to as 

pseudospins, since they result in the single particle wavefunctions having a spinor form pp. 

In these proceedings we examine a monolayer of graphene in the presence of a strong 
perpendicular magnetic field of magnitude B. The Landau level (LL) energies for graphene 
are, e n = sign(n)/iVF y / 2|n|/^s, where vp ~ 10 6 ms _1 is the Fermi velocity, £g = \Jhc/eB is 
the magnetic length and sign(O) = 0. There is also the usual LL degeneracy resulting from 
the different possible guiding centres for electronic cyclotron orbits. This is described by the 
oscillator quantum number, m, in the symmetric gauge. We consider magnetoplasmons (MPs), 
which are essentially excitons created from a ground state described by the filling factor v = 2, 
where all states with LL indices n < are occupied and all other states are empty. An exciton 



is formed when one of the electrons in this ground state is excited to a state in a higher lying 
LL and becomes bound to the hole it leaves behind. We focus on excitons with the lowest 
possible excitation energies, namely those with the hole in the n = LL and the electron in the 
n = 1 LL. Taking into account the spin/pseudospin states of the electron and hole, this yields 
16 resonant excitations. These are collective excitations on account of the LL degeneracy and 
also due to the fact that excitons with different spin/pseudospin characters may be mixed by 
the two-body Coulomb interaction. In the following, we study the SU(4) symmetry in graphene, 
which arises from the two spin and two pseudospin degrees of freedom, denoted by f, ! and ff, 4 
respectively. This leads to a direct analogy between MPs in graphene and mesons composed 
of either first or second generation quarks. We then review the dispersion relation for MPs in 
pristine graphene and finally examine the states which become localised in the presence of a 
single Coulomb impurity. In both cases we observe the degeneracies predicted by the SU(4) 
symmetry and Young diagram techniques. 



2. SU(4) symmetry in graphene 




Figure 1. SU(4) [15]-plet describ- 
ing the excitonic states in terms 
of flavour isospin (S 2 ), hypercharge 
(Y) and charm (c). 
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Figure 2. Dispersion relations for 

excitons with the hole in n = LL 

and electron in n = 1 LL. V is the 
quasimomentum in units of h/£s- 



An electron (or hole) in graphene in a strong perpendicular magnetic field is described 
by four quantum numbers: the LL index n; the oscillator quantum number m; the spin 
projection s z and the pseudospin projection t z . We denote the fermionic creation operator 
for an electron by c nm T Z s z and that for a hole by d nm T Z s z = c nmTzSz . The first two orbital 
quantum numbers are associated with the cyclotron orbits. The last two, when considered 
together, give four possible spin/pseudospin "flavours", leading to the SU(4) symmetry. We 
identify these with the four quark flavours up (u), down (d), charm (c) and strange (s) according 
to {JJ4-, 1H, JJ-t> fir} = {d,u,s,c}. The generators of SU(4) can be constructed using bilinear 
combinations of the fermionic creation and annihilation operators conserving the number of 
electrons and holes. We have 

A? = ^2 C nmi C nmj ~ ^ d) nm -d nm i, (1) 
nm nm 



where i, j G {d, u, s, c}. The generators satisfy the commutation relations 



[Cij, Cki] — 8jkCu — duCkj, 



(2) 



so are closed under commutation, confirming that they form a Lie algebra. The operators 
associated with spin and pseudospin may of course be expressed in terms of the generators. For 

example, S z = \ (c ss + C cc - C uu - C dd ^j and S + = S x + iS y = h (c cu + C sd ^j . 

In particle physics, mesons, which are a bound quark-antiquark pair, are classified using 
Young diagrams [6]. We use the same techniques to determine the multiplet structure for 
electron-hole (e-h) complexes in graphene. Young diagrams describe SU(n) multiplets and 
consist of n or fewer left-justified rows of boxes, such that each row is as least as long as 
the row beneath it. Roughly speaking, each box represents a particle. Each SU(n) multiplet is 
uniquely labelled by a list of n — 1 non-negative integers, (ax, a 2 , . . . , a n _i). The integer in the 
i th position corresponds to the number of boxes in the i th row minus the number of boxes in 



the (i + 1) row i.e. the overhang of the i row. By the definition of a multiplet label, adding 
columns of length n to the left of a Young diagram is redundant. Young diagrams can be used to 
calculate the multiplicity, [iV(n)], of an SU(n) multiplet. For the n = 4 case of current interest, 



[JV(4)] 
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(3) 



There are two fundamental representations of SU(n): [n], with Young diagram and label 
(1,0, . . . 0) and [n] with Young diagram composed of n — 1 vertically stacked boxes and label 
(0. . . ,0, 1). The latter represents an antiparticle; this makes sense if we identify n vertically 
stacked boxes with the vacuum. We shall also use the symbol, 1 to denote the vacuum. For 
the current purpose, the hole in a MP behaves as an antiparticle. Thus the electron has Young 



diagram and the hole Young diagram 



To determine the SU(4) multiplet structure for 



the MPs, we must combine these two diagrams. The rules for combining Young diagrams [6] 
yield 



1. 



(4) 



Eq. @ is equivalent to (1,0,0) ® (0,0,1) = (1,0,1) (0,0,0) or [4] ® [4] = [15] [1]. Now 
the generators, Cij, commute with the interacting Hamiltonian describing the electron system 
up to some small symmetry breaking terms [Jj- In addition, the application of a generator to 
a state in a particular multiplet always returns a state in that same multiplet with the same 
Coulomb interaction energy. This ensures that states within a given multiplet are degenerate 
in their energies. Thus from this simple calculation, we have shown that the MPs in graphene 
will either have degeneracy 1 or 15. The [15]-plet is shown in Fig. [I] and is analogous to the 
meson multiplet of particle physics. The states are plotted as a function of: the flavour isospin 

projection t z = \ {c uu - C d d}, the hypercharge Y = \ (c uu + C dd + C c ^j - §(7^ and the 

charm c = C cc . One may deduce the symmetry of states under exchange of quantum numbers 
for different particles from their Young diagrams [6]. 



3. Dispersion relation for magnetoplasmons in pristine graphene 

As a precursor to our results for localised MPs in dirty graphene, we discuss the dispersion 
relation for extended MPs in a clean system. The calculation is explained in Ref. [8] in the 
same spirit as that for the 2DEG established by Kallin and Halperin [§], but with the additional 




Figure 3. MPs from the ground state v = 2 localised on (a) a donor impurity and (b) an 
acceptor impurity. The hatched area of width 0.75-Eo represents the continuum of extended 
MPs. Energies are given relative to the lower continuum edge in units of Eq as a function of the 
generalised angular momentum projection, J z . 



complication of the valley pseudospin degree of freedom. The centre of mass quasimomentum 
for the exciton, P, can be shown to be a good quantum number. One can think of this as being 
due to the repulsive Lorentz forces and attractive Coulomb forces, which act on the electron and 
hole, cancelling out when they have a certain separation, enabling the exciton's centre of mass 
to move in a straight line. The Hamiltonian (in the symmetric gauge) for MPs is |10j : 

Af,m Af,m 
+ ^2 ^2 ^^^miV^m^A^m^jv^m^^maCMmi, (5) 

where we define the collective index N = {nT z s z }. The first two terms give the single 
particle contribution. The e n are the LL energies renormalised by a self energy correction 
due to exchange interactions. The last term describes the e-h interaction; it has direct and 
exchange components. This Hamiltonian mixes the four excitons for which there are no spin or 
pseudospin flips and leaves the remaining twelve unmixed. The Coulomb interaction is treated 
exactly in the span of given LLs and the dispersion relations are found by diagonalising the 
Coulomb interaction Hamiltonian in the basis of two-particle electron-hole states compatible 
with magnetic translations (having definite quasimomentum P). The Coulomb interaction is 
treated as a perturbation and the dispersion relation found by diagonalising the Coulomb 
interaction matrix in the basis of these non-interacting two-particle wavefunctions. For the 
fully filled LL considered here, there are only two possible dispersion relations. The four mixed 
excitons give a non-degenerate dispersion branch [H] 



(6) 



and a triply degenerate dispersion branch 



E 2 = Ex + E^Ve- — , 



(7) 



where V = Pis/fr is the dimensionless centre of mass quasimomentum, P, for the exciton and 
Iq and I\ are modified Bessel functions. The dispersion, E2, is the same as for tthe 12 unmixed 
excitons, making the total degeneracy 15, as predicted in the previous section. These dispersion 



relations are plotted in Fig. 
Coulomb interaction energy, 



Energies are given in units of E$ — y 2 elB . 
where e is an effective dielectric constant. 



the characteristic 



4. Localised magnetoplasmons in dirty graphene 

We now consider how MPs may become localised in the presence of a single Coulomb impurity 
with potential, V(r) = where Z is the impurity charge in units of e [10]. The centre of 
mass quasimomentum is no longer well defined for a MP. However, the axial symmetry remains 
and the orbital angular momentum projection, defined in terms of the electron and hole orbital 
quantum numbers as, J z = \n e \ — \nh\ —m e + m,h, is a good quantum number. In the presence of 
an impurity, localised MPs break away from the continuum of extended MPs as discrete states. 
These localised MPs are shown for a donor [Z = 1) and acceptor [Z = —1) impurity in Fig. 
[3j plotted as a function of J z . The yellow shaded area represents the continuum; it has width 
0.75-Eb, as can be seen in Fig. [2j Note that for both impurity cases, there are non-degenerate 
branches and branches with 15-fold degeneracy. The predictions based on SU(4) symmetry hold 
here, because the Coulomb impurity, being long ranged, does not introduce an inequivalency 
between the valleys. This would not be the case for a short-ranged impurity located on one of 
the carbon sites |llj . 



5. Conclusions 

In conclusion, we have studied neutral collective excitations of graphene in the presence of a 
strong perpendicular magnetic field. Electrons and holes in graphene can exist in four possible 
spin/pseudospin states, so that electron-hole complexes may be likened to mesons composed of 
quarks with four possible flavours. This analogy enabled us to use techniques established for 
treating mesons to obtain the multiplet structure of the neutral collective excitations for both 
extended excitations and those localised on a Coulomb impurity. An extension of this idea to 
charged collective excitations, which can be thought of as three particle complexes, can be found 
in Ref. [12]. These are analogous to baryons formed of three quarks. 
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